EXPONENTIAL APPROXIMATIONS FOR THE 
PRIMITIVE EQUATIONS OF THE OCEAN 

R. TEMAM AND D. WIROSOETISNO 



Abstract. We show that in the limit of small Rossby number e, the primitive 
equations of the ocean (OPEs) can be approximated by "higher-order quasi- 
geostrophic equations" up to an exponential accuracy in e. This approximation 
assumes well-prepared initial data and is valid for a timescale of order one 
(independent of e). Our construction uses Gevrey regularity of the OPEs and 
a classical method to bound errors in higher-order perturbation theory. 



1. Introduction 

We consider the primitive equations for the ocean (henceforth OPEs), scaled as 



in 



d t v + i [v x + Vp] + u ■ Vv = t iA 3 v + f v , 



(1.1) dtp-jw + U'Vp = [i,A 3 p + f p , 
V • u = V • v + w z = 0, 

P = -Pz- 

Here u = (u, v, w) is the three-dimensional fluid velocity, with v — (u, v) its hori- 
zontal component and v — (—v, u); p is the pressure; p is the perturbation density 
(not including the mean stable stratification which figures into the e in the equation 
for p). We write V2 := (d x ,d y ) and V3 := (d x ,d y ,d z ); when no ambiguity may 
arise, we simply write V. Similarly, we write A 2 := d% + dy and A 3 :— d^+dy+d^. 
The parameter e is related to the Rossby and Froude numbers; in this article we 
shall be concerned with the limit e —> 0, and for convenience we assume that e < 1 
(further restrictions on e will be stated below) . In general the viscosity coefficients 
for v and p are different; we have set them all to p for clarity of presentation (the 
general case does not introduce any more essential difficulty). The forcings f v and 
f p are assumed to be independent of time. 

We work in three spatial dimensions, x = (x, y, z) 6 [0, Li] x [0, L2} x [— L3/2, L3/2] 
=: M, with periodic boundary conditions assumed. Following common practice in 
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numerical simulations of stratified turbulence (see, e.g., 0), the dependent vari- 
ables are assumed to have the following symmetries: 

v(x,y,-z) = v(x,y,z), p(x,y,-z) = p(x,y,z), 

w(x,y, -z) = -w(x,y,z), p(x,y,-z) = -p(x,y,z); 

we say that v and p are even in z, while w and p are odd in z. If in addition 
f v is even and / is odd in z, it can be verified that this symmetry is preserved 
by the OPEs (D) , that is, if it holds at t = 0, it continues to hold for t > 0. 
Since w and p are periodic in z, we have w(x,y, — L3/2) = w(x,y, L3/2) = and 
y, — L3/2) = p(x,y, L3/2) = 0; similarly, w z = 0, v z — and p z = on 
z = 0, ±£3/2 if they are sufficiently smooth (as will be assumed below). One may 
consider the symmetry conditions Ijf ,2|l as a way to impose the boundary conditions 
to = 0, p = 0, u z = 0, v z = and p z ~ on both 2 = and z = L3/2, in the effective 
domain [0,Li] x [0, L2] x [0, L3/2]. All variables and the forcing are assumed to 
have zero mean over M; the symmetry conditions above ensure that this also holds 
for their products that appear below. 

It is known that, given sufficiently regular initial data, the OPEs have a unique 
strong solution for all time 0],jZj; the results of ^3] then imply that a very regular 
solution, belonging to a Gevrey space defined below, exists for all time. The exis- 
tence of finite-dimensional global attractors in various spaces (see QHDj an d thus 
uniform boundedness of the solution, have been shown in |H] ; see also the result in 
|12j . One can thus regard the OPEs as equivalent to a finite-dimensional system. 
This fact is an important ingredient for our main result. 

As noted earlier, we are concerned with the limit e —> 0. A flow is said to be in 
geostrophic balance if the 0(1/ e) terms in ((1 . lh ,^ vanish, i.e. if 

(1.3) v = V ± p 

where V :— (—d y ,d x ). Note that this implies 'V ± .v = 0, which with the fact 
that w is odd in z in turn implies that w = 0, so the 0(1/ e) term in l|f ,1b ) also 
vanishes. Assuming that (|I.3() is satisfied initially, it is well known that the OPEs 
can be approximated by the simpler quasi-geostrophic equation (QGE), 

(1.4) d t qZ + vZ-VqZ = pA 3 qZ + f q , 

which only involves a single variable q g ; here / := ■ /„ — d z f p . The original 
variables are recovered using 

(1.5) v e = V ± Az 1 q s and p s = -^A^V, 

Here A^ 1 is uniquely defined to have zero mean over M. The quantities V s and 
p g derived from the potential vorticity q g are said to be geostrophic. The fact that 
d z w + V • v = and (fTBT ^ imply that w s = 0. 

An approximation result was obtained in |5j for a closely related system where it 
was shown that the QGE (|I.4|I with p = and f q = approximates (the unforced, 
inviscid and Boussinesq analogue of) the OPE (|f .I|) with an error of order e for 
t G [0, T], viz., 

(1.6) \v(t) - (t)| ff4(M) + \w(t)\ HHM) + \p(t) + d z A^ q e(t)\ HHM) < e 

assuming that the left-hand side less than e at t = (and given sufficient regular- 
ity). We shall not attempt to follow their approach in this article as it makes no 
qualitative difference for our exponential-order result. 
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Regarding this result as a first-order approximation in the parameter e, a natural 
question is whether one can obtain higher-order approximations. Put differently, 
we would like to know how fast, as £ — > 0, the solution of the PE converges to 
the solution of a simpler system analogous to the QGE l|1.4|l . Our main purpose 
in this article is to show that convergence of any algebraic order is possible for e 
sufficiently small, resulting in an error estimate which is exponentially small in s. 

In the geophysical parlance, the dynamics of a fluid flow is said to be balanced if 
the solution stays near a subspace ("balance manifold") for some time. The order 
of the balance dynamics measures how close this approximation is in terms of e. A 
result such as ()1.6|l implies that the quasi-geostrophic relations H1.5(l define a balance 
dynamics of order e. Our result (Theorem ^ below) thus implies that a balanced 
dynamics of exponential order exists for a timescale of order one (i.e. independent 
of the Rossby number). Note, however, that our result is actually stronger than 
the geophysical definition of balance, since it gives a pointwise approximation to 
an exponential order. 

Also of geophysical interest is whether the solution (when appropriately ini- 
tialised) stays close to a balance manifold over longer timescales. Using careful 
estimates (of a different type than those used here), it may be possible to prove 
that the solution remains "exponentially balanced" for timescales of order 1/e (al- 
though it is clear that one cannot expect pointwise accuracy over this timescale); 
we plan to report on this in a future work [lfij . 

We begin with some notations. We write the horizontal velocity v as 

(1.7) v(x,y,z) ^^(zJ + V^ + Vx 

where v(z) is the mean vertical shear, i.e. the (x, y)-average of v. We introduce the 
linearised potential vorticity (which has appeared in the QGE above) 

(1.8) q := V 1 - -v-p z . 

The streamfunction ip and velocity potential x are defined as follows. First, let 



where, here and henceforth, A 2 is uniquely defined to have zero mean on each 
(x,y)-plane. Next, let 



where A3 is uniquely defined to have zero average over M. This completes the 
definition of (q, v, x, <fi) in terms of (v, p). 

It can be verified that these definitions imply the following relations, 

(1.12) A 2 V> = V- L v, p = ijj-ct> and A 2 x = V • v. 

The second equation is useful to compute p = —d z p from q and <f>; v can be com- 
puted from (q, v, x, 4 1 ) using (|1.7(l and Hl.llfl . Using the incompressibility condition 
V • v + w z = and the fact that w(x, y, 0) = 0, we also have 



(1.9) X -=^V-v 



Jo 

And finally, let 




(1.11) A3 1 + 
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The primitive equations can now be written in the following form, 
t q + V x • (u • Vv) -8 z {u- Vp) = pA 3 q + f q , 

d t v z + -Vg + d z (wv z ) = fiv zzz + d z f v> 

e 

(1-14) d t A 3X -^A 3( p + A 3 A^V-(u-Vv) = pA 2 3X + A 3 f x , 

d t <t>zz + ^A 3 x + d^A^V 1 - ■ (u ■ Vv) + P z d z {u ■ Vp) 

= pA 3 <p zz + d zz f c j > , 

where overbar denotes average over (x, y) and P z is a projection which removes 
the (x, y)-average of a function, P z p := p — p, etc. The forcing terms in (|1.14|) are 
related to those in by 

f q = V x -f v -d z f p , f x = A^V-f v , 

L=Z, dzzh = ^'dzzV 1 - ■ f v + P z d z f p . 

It is easily verified that the 0(l/e) terms in Ijl.l4|l are antisymmetric. If the fast 
variables (v,X> 4 1 ) are a H zero — a condition that is equivalent to the geostrophic 
balance relation — they can only grow through the nonlinear terms involving 
the slow variable q and the nonlinear terms on the right-hand side. Further analysis 
(cf. [3],[J§]) shows that they remain small over a timescale of order one. Our main 
result [Theorem ^ states that, for certain initial data (v, \ an d 4> given in terms 
of q), this system can be approximated up to an error exponentially small in e by 
the solution of an equation for a finite-dimensional slow variable q. 

We note that the form (|1.14|) will only be used for (i) the construction of the 
approximation v*(q; e) and p*(q; e), and (ii) the integration of the finite-dimensional 
variable q. The approximation error estimates are in terms of (v,p) and (v*,p*), 
so we will not need a priori estimates for (|1.14(l . 

2. Gevrey Regularity and Statement of the Main Result 

Let W = (u, v, p). The periodic boundary conditions allow us to write functions 
(assuming sufficient regularity) in Fourier series, 

(2.1) W(x,t)= J2 W k (t)e ik ' x , 

where k' := (k[, k' 2 , k' 3 ) with k[ := 2ir ki/L{. For a > fixed, the Gevrey norm || • Ho- 
is defined as 

(2.2) \\W\\l := e 2ff|fc| |VF fc | 2 . 

fcez 3 

It follows from this that if HWHo- < oo, \Wk\ decays exponentially fast in |fc|; this 
fact will play an important role below. 

The following result on the OPEs is proved in (this has been extended to all 
time in |12j . but the version here is more useful for our present need): 

Lemma 1. Let Wq = (uq,vq, po) be such that IV3VF0I = C\fl < 00 and sup- 
pose that \V 3 f v \ + IV3/J = Cf < 00. Then for any C\ > Ci,o there exists a 
Ti(Ci )0 ,C/, d;(x,fi,M) > such that for t e [0,Ti] ; 

(2.3) \V 3 W(t)\<C 1 . 
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Now let a > be fixed, and let Wo be such that \\A 3 Wo\\<j = C^o < oo; in addition, 
suppose that || V3/„||cr + || ~V 3 f p \\a = C't < oo. Then for any C a > C a ,o, there exists 
a T a (C<jfl, C'f, C a \ a, /i, M) > such that for t £ [0, T a ], 

(2.4) \\^W{t)\\ a <C a . 

We stress that T\ and T a are independent of e. (Here and henceforth, all constants 
are understood to be positive.) 

Given a fixed k > 0, we define a finite-dimensional truncation of W as 

(2.5) W<(x,t) := (P'WOOM) := £ TU fe (i)e ifc '*, 

|fc|<K 

and W > :— W — W < . It is clear from this definition that the projection P < is 
orthogonal. For the low modes we have the "reverse Poincare inequality" , 

(2.6) \V 3 W<\ S <cn\W<\ s . 

Here and throughout this article, || ■ || CT denotes the Gevrey norm (E3 and | • | s 
denotes the usual Sobolev H s norm; when no confusion may arise, we will often 
write simply | • | for the L 2 norm | • | . If ||W||<7 < oo, the exponential decay of the 
Fourier coefficients Wk allows us to bound the high modes W > in L 2 as 

(2.7) \W>\o<C^-™\\W\\ a 
or, more generally in H s , 

(2.8) \W > \,<C.K'&-'" t \\W\\ a . 
We can now state our main result: 

Theorem 1. Let a 6 (0, 1) be fixed. Suppose that the initial data Wo — (uo, i»o, Po) 
satisfies 

(2.9) HViWbL < oo, 
and the forcing satisfies 

(2.10) ||Vi/ B || <r +||V|/ p || cr < o. 

Then there exists an £o(|[^^3'Wo||er) 11^3/llo-i a i Mj > 0, and a q (which is finite 
dimensional) evolving according to 

(2.11) dtq+P" [V x ■ (u* ■ W) - d z (u* ■ Vp*)] = pA 3 q + f< , 

where u* = u*(q;e) and p* = p*(q;e) are constructed in the proof below, such that 
if e € (0,£o) an d the initial data satisfies 

(2.12) K - »*(g(0);e)|g + \ P o ~ P*(g(0);e)| 2 < C id exp(-2 ( r/ £ 1 / 4 ), 

for some constant C;d {independent of e and Wo) then a similar estimate, 

(2.13) \v(t)-v*(q(t);e)\ 2 o + \p(t)-p*(q(t);e)\ 2 < 4C id exp(-2a/e 1 / 4 ), 
holds fort e [0,T„] where T* = T*(\\VlW \\ a , \\ V 3 3 f\\ ff ; a, fi,M). 

Remarks. 

1. The v* and p* are the "exponential order" analogues of the (leading-order) 
quasi-geostrophic v g and p s . Unlike the latter, which is infinite-dimensional, our 
higher-order approximations are based on a finite-dimensional variable and depend 
on Gevrey regularity of the parent system (|1.1|) • 
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2. The sharpness of our bound (|2.13(1 is unclear, but an explicit example con- 
structed in ^2] for a closely related problem suggests that one cannot do better 
than exp(— c/e). 

3. Although in this article we only treat the primitive equations for the ocean 
Ijl.lfl . the method described here should be applicable to many PDEs with a small 
parameter for which a Gevrey regularity result can be proved. 

3. Proof of Theorem [I] 

Our approach is inspired by which is an averaging- type result on single- 
frequency infinite-dimensional systems, and by [SJ App. B], which considers a finite- 
dimensional problem. Our present problem is infinite-dimensional with an infinite 
number of frequencies, so small denominator problems would appear in an extension 
of Jl] to this case. However, by considering the singular perturbation problem 
(hence the requirement for special initial data) in this article, we can avoid the 
small denominator problem. 

The plan of the proof is as follows: 

* We will build up our nth-order approximation W n as: 

(3.1) W = W< + W> 

(3.2) W<=W n + W 

The high modes W > are exponentially small thanks to the Gevrey regularity of 
W; (|2.8I) then implies that they are exponentially small in any Sobolev norm. The 
low modes W < (which are finite-dimensional) can then be approximated by W n , 
which is constructed such that W is exponentially small in e. The total error 
W — W n — W > + W is then also exponentially small. Like W < , the approximate 
solution W n = (v n ,p n ) is finite dimensional, 

(3.3) W n (x,t)= J2 Wftfje*-. 

|fc|<K 

* It is "slaved" to the slow variable q in the following manner: 

(3.4) v n — V n + V-V" + VX" with A 3 V>™ = q + d zz <S> n 

(3.5) p"-A 3 - 1 [A 2 9 z $"-a 2 g], 

where V(q;e), X n (q;e) and $ n (q;e) are functions to be computed below. For a 
motivation of this "slaving ansatz" , we refer the reader to |20| . Even when work- 
ing with finite-dimensional systems, we shall keep the convenient PDE notation, 
remembering that in this case we have the reverse Poincare inequality (|2.6(l . We 
note for future reference 

(3.6) w n = -f A 2 X n (g;e)dz'. 

* The (truncated) potential vorticity q evolves according to our "higher-order quasi- 
geostrophic equations" (QGE n ), 

(3.7) d t q + P" [V x • (u n ■ Vv n ) - d z (u n ■ Vp n )] = nA 3 q + f< , 

where u n — (v n ,w n ) and p n are determined in terms of (q, V™, X n , $ n ) by (|3.4|l - 
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* In turn, (V" ,X n ,$ n ) are computed using the iteration (|3.11j) - (|3.13j) . If this 
iteration were convergent, say to (V 00 , X°° , $°°), and if we let q be the solution of 

(3.8) d t q+ P < [V x • (tt 00 • Vu°°) - d z (u°° ■ Vp°°)] = ^A 3 q + f< , 

where it 00 and p°° are determined from (V 00 , X°°, $°°) as before, then (v 00 ^ 00 ) 
would be an exact solution of the OPEs. As is often the case in this type of problems 
(see ^U] and ^7] for closely related problems), this iteration is asymptotic rather 
than convergent, so we have to end the iteration at some n = n*. It is shown 
below that for n = 0, • • • , n* = \r\je x ^\ : (i) the slaved variables (V™, X™, $") are 
bounded [see 1|3. 39(1 below], and (ii) the QGE™ has an error of order e n locally in 
time [see ()3.56|) below]. 

* We then show that with v n = V n (q;e), X n = X n (g;e) and 0" = $"(g;e) the 
solution of i|3.7[l is bounded for < t < T with T independent of e. (One can in 
fact show that the solution is bounded for all time, but we shall not do so here.) 

* Finally, we obtain bounds for W using the usual Gronwall-type argument. This 
proves the main theorem since W > is exponentially small by Lemma ^ 

Following a common practice, we write c for a generic constant which may not 
be the same each time it appears; the more important constants are numbered: c x , 
c 2 , etc. Unless otherwise indicated, these constants (assumed positive) may depend 
on s, a and domain size (Li), but not on s, k, n, p or the initial data. 

3.1. Construction of QGE". Let V° = 0, X° = and $° = 0. From fi ty -fi ty 
we then have 

(3.9) v° = V^A^q and p° = -d z A 3 1 q. 

Putting these in (|3~7jl and using the fact that V ± -(v°-'Vv )-d z (v -Vp ) = v°-Vq, 
we find 

(3.10) d t q + P" [v° ■ Vq\ = pA 3 q + f< , 

which is the quasi-geostrophic equation (|1.4|l - (|1.5|l for the truncated variable q. 
This is our zeroth-order (truncated) quasi-geostrophic model, QGE . 

For n > 1, we define V (q;e), X n (q;e) and <fr n (q;e) by the following iteration 

[cf. irrat ] 

(3.H) - \{VTY = P < {(DV^ (q, V n , X", $") + d z (^™) } - pd zz V n z - dj v , 

1a 3 $™ +1 = p < {(DA 3 x n )g < {q, v n ,x n ,$ n ) 

(3.12) + A 3 A 2 - 1 V • (u* • W 1 )} - pA 2 3 X n - A 3 f<, 
--A 3 X n+1 = P < {(DA 3 $") g K (q, V n X n , $") 

(3.13) + d zz A^V x ■ (u n ■ Vv n ) + P z d z (u n ■ V P n )} - pA 3 ^ z - d zz f< 7 
where D denotes derivative (differential) with respect to q and 

(3.14) Q < (q, V n , X", $") = P" [fiA 3 q + f q - V x • (u n ■ Vv n ) + d z (u n ■ Vp n )] 

= :g{q,V 7 \X n ,^) + ^A 3 q + f<. 
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We recall that (u n ,p n ) are given in terms of (V™,X™,$ n ) by (23)) -(EH). This 
construction was presented, in general and formally, in |18) . which also proposed 
the use of a series expansion. 

Let s > 3/2 be fixed (in the theorem we use s = 2, but we keep the general s 
in the proof), and consider q fixed for now. We start by estimating the nonlinear 
terms in Q < (q, V n ,X n , $") =: Q n using (fT"TU|) and ^Z$-&f$: For d z (u n ■ Vp n ), 
we compute 

KI,<c(|y n | s + |v 2 v/ l | s + |v 2X n | s ) 

< c(\V n \ s + \q\ s -! + |V 3 $"| S + |V 3 X"| S ), 

. ^ \v% < c(\V n z \ s + \V 2 ^\ S + \V 2X %) 

(3.16) __ 

<c(\V n z \ s + \q\ s + \A 3 <!> n \ s + \A 3 X n \ s ), 

(3.17) K| 5 <cK| s <c|A 3 X"| s . 
From these we have 

c\u n \ a < \V\ + |gU_! + |A 3 X% + |V 3 $"| S 
c|<| s <|y ? ;| s + |g| s + |A 3 X"| s + |A 3 $"| s . 
Putting these together with 

(3.19) c\V 3 p n \ s <\q\ s + \A 3 $ n \ s and c\V 3 p n z \ s < \q\ s+1 + \A 3 <P n \ s+1 
gives us 

|< ■ Vp n \ s < c (\q\ s + \V% + |A 3 $"| S + |A 3 X"| S ) (\q\ s + |A 3 $"|.) 

|«" ■ Vp«| s < c (|qU_i + |V n | s + |V 3 $ n | s + |A 3 X"| s )(|g| s+1 + |A 3 $"| S+1 ) , 
and, upon the use of H2.6(l . 

(3.20) \8 z (u n z ■ Vp n )\s < en (\q\ s + \V% + |A 3 $"| S + |A 3 X"| S ) 2 . 
For the other nonlinear term in Q n , 

(3.21) V 1 - ■ (u n ■ Vv n ) = u n ■ VA 2 f + {V ± u n ) : V 3 {V n + V ± t/j n + VX"), 

we need in addition 

|V 2 it"| s < c(\A 2 ip n \ s + |A 2 X n | s + \V 2 w n \ s ) 
(3 ' 22) <c(|g| s + |A 3 X"| s+1 + |A 3 $"| s ), 

\V 3 v n \ s <c(\V n z \ s + \A 3 r\s + |A 3 X"| S ) 

<c(\q\ s + \V n z \ s + \A 3 X n \ s + \A 3 ^\ s ), 

and 

(3.24) c | V 3 A 2 <T| S < |g| s+1 + |A 3 $"| S+1 . 
Putting these together and using l|2.6|) again, we find 

(3.25) IV 1 - • (u n ■ Vv n )\ s < ck (\q\ s + \V n z \ s + |A 3 $"| S + |A 3 X"| S ) 2 . 
The linear terms are easily bounded, giving us 

(3.26) |<T| S < Cl « (\q\ s + \V n z \ s + |A 3 $% + |A 3 X"| S ) 2 + [l K 2 \q\ s + |/<| s . 
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) 2 

s > 



Going next to the nonlinear terms in (|3.II|) - (|3.I3|1 . we compute 

(3.27) \W^\ 8 <c\w n \ s \V x ^ + VX% < CK |A 3 X"| s (|g| s + |A 3 $"| s + |A 3 X 

The following estimates follow from our computation above: 

IA3AJ 1 V • (u n ■ Vv n )\ s < ck* (\q\ s + \V n z \ s + |A 3 $"| S + |A 3 X» 
(3.281 _ „ 

Id^A^V 1 - ■ (u n ■ Vv n )\ s < ck 3 (\q\ s + \V n z \ s + |A 3 $"| S + |A 3 X"| S ) 2 . 

Using the relations Ijl.l5|l . the forcing terms are bounded as 

|/,|.<c(|/.Ui + |/ p |.+i) |A 3 / X | S < c\f v \ s+3 

\dzfv\s < 1/oL+l \9zzU\ s < c(\f v \ a+3 + |/p| s+3 )- 

Moving on to terms such as D(A 3 X n )Q n in l|3.11[l - i|3.13|) . we need a few defi- 
nitions: Let D v (q) be the complex ^-neighbourhood of q in P iP(M). It can be 
defined by Fourier series as 

(3.30) D v (q) = [q' : q'(x) = £ q'J kx with ]T |fc| 2s \q' k q k \ 2 < V 2 } . 

\k\<K \k\<K 

We note that since q{x) is real, the Fourier coefficients must satisfy q-k = q\ 
where (here and only here) bar denotes complex conjugate; however, the Fourier 
coefficients q' k in l|3.30|l need not satisfy this constraint. Let 6 > be fixed. For 
any function g of q, let 

(3-31) \g\ s . n := sup \g(q)\ a , 

Dn-ns(q) 

which is meaningful for n = 0, • • • , \jl/6\ =: when Z) I) _„ ( 5(q ; ) is non-empty. 
Let us also denote 

(3.32) \U n \ s , n := \V n z \ s , n + |A 3 $"| s; „ + |A 3 X"| S; „. 

To bound D(V")G n , D(A 3 X n )G n and D(A 3 $™)£™, we use Cauchy's integral for- 
mula (cf. e.g., [§1): For <p : c\D v (q) — » C analytic with c\D v (q) C C' and S > 0, 
one can bound \ip'\ in c\D n -s(q) by \ip\ in clD^q), 

(3.33) |y>' ■ z| clDrj _, te) < ^Mcl £>„(§) Mc«- 

Now V" z (9;e), A 3 X"(g;e) and A 3 $™(g;e) are analytic functions of q, being poly- 
nomials in q k . Using this, we have 

|(DA 3 $")S"| S! „ +1 < -AA^ n \ s Jff\ n 
(3-34) J 

< -|A 3 $"| 5 , n { M « 2 |g| s + !/,<!. + ™ 2 (|g| a + |c/"| s; „) }• 

We have now bounded every term in l|3.13|l : putting together (|3.34() . I|3.29|l . 
and ipOH)l . we find 

i|A 3 X" +1 | s;n+1 < /iK 2 |A 3 $"| S .„ + |a„/<| a + ck 3 (\q\ s + \U n \ s , n ) 2 
(3.35) £ 

+ ^\A 3 ^ n \ s , n {^ 2 \q\ s + \f<\ 8 + ck 2 {\q\ s + \U n \ sin ) 2 }. 
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Similarly, we find for (|3.12|1 . 

-\A 3 <P n+1 \ s ., l+1 < /in 2 |A 3 X"| S; „ + |A 3 /<| S + ck 3 (\q\ s + \U n \ s ., n y 

(3.36) e 

+ - 5 \^ n \ s , n {^ 2 \q\ s + \f<\ s + CK 2 (\q\ s + \U n \ s , n f}, 

and for (|3.11|) . 

-\V r : + \ n+1 < [IK 2 \V n z \ s , n + \dJ v < \ s + CK 3 {\q\ s + \U n \ s ., n ) 2 

(3.37) £ 

+ C -\A^ n \ s , n {fl K 2 \q\ s + \f<\. + CK 2 (\q\ s + \U n \ s , n ) 2 }. 
Putting these estimates together with (|3.29|) . we find 

-\U n+1 \ svn+1 < yLK 2 \U n \ s;n + ll/H + C 2 K 3 (\q\ s + \U n \ s - n ) 2 

(3.38) £ 

+ | \U n \ s . n {nn 2 \q\ s + \f<\. + c 4 « 2 (|g| a + \U n \ sin ) 2 }, 

where ||/|| := |(/<, /<)| 5+3 . 

Now suppose that at iteration n, with c x = 8c 2 we have 

(3.39) \UXn<Ci£ 1/4 (\€ + \\f\\)- 

Since V° = X° = $° = 0, this trivally holds for n — 0. Taking S = e 1 / 4 and 
K = e -1 / 4 , and substituting this into (|3.38|l . we find that 

(3.40) |C/" +1 | s;Ii+1 <c l£ 1 / 4 (|gt + n/ll). 
provided that e satisfies the following constraints 

e <min(l/(16/i 2 ),l/ C 4 ! ( Cl /4) 4 / 3 , 1 / ( 1 2 M c 3 1 ^ | s ) 4 , 

( 3 - 41 ) 

l/(12c 3 ||/||) 4 / 3 , l/(l2 C3 c 4 [2|g| s + ||/||] 2 ) 4 ). 

By induction (|3.39(l can be seen to hold for n = 1, • • • , n*. For future reference, we 
note that this restriction on e and l|3.39|l imply that 

(3.42) \q\ a + \U n \.<c 3 (\q\ a + 



3.2. Local Error Bounds. Moving on to estimating the approximation error lo- 
cally in time, we define 

(3.43) n% : = (Dv n z )g n + -v n z + d z (W^) - fid zz v: - dj< 

£ 

= i[v:-v: +i ] 

(3.44) Tl n x := (DA 3 X")£ n - -A 3 $" + A 3 A^V • (u n ■ Vv n ) - a*A 2 X" - A 3 /< 

= J[A 3 $" +1 - A 3 $"] 

(3.45) Tq, := {V<S> n zz )G n + Ja 3 X" + d zz A^V x ■ (u n ■ Vv n ) 



+ P z d z (u n ■ Vp n ) - M A 3 $S, - d zz f< 
- [A 3 X" - A 3 X' i+1 ] . 



EXPONENTIAL APPROXIMATIONS FOR PRIMITIVE EQUATIONS 11 

Since V° = 0, X° = and $° = 0, we have 

n° x = A 3 A^V-(v°-Vv°)-A 3 f<, 
K° = d zz A 2 l V x ■ (v° ■ Vv°) + P z d z (v a ■ Vp°) - d zz f<, 
which leads to the estimate, denoting \TZ n \ s := |ftg| s + [R*\ s + \W£\ S , 

(3.46) |^ | s <C £ - 3 / 4 (k1 s + ||/||) 2 . 

For n > 1, let us denote 5u™ := — v n , Sp n := p n+1 — p n , etc., and consider 

= (DV n z +1 )G n+1 + -VT 1 + d z {w^v n z +1 ) - pd zz V n z +l - dj< 

(3-47) _ (DV' z +1 )5g n - e{DKl)g n 



+ d z {8w n vf) + d z (w n+1 8vV;) + epd zz K™ , 

For the second equality, we have used the second equalities in H3.43|l - (|3.45|) to write 
5v n in terms of 1Z%, etc. To estimate these 5-quantities, we compute 

A 3 Sp n = A 2 [d z S<P n ] 
(3 ' 48) => \V 3 Sp n \ s <c\A 3 S^\ s <ec\1l-\ s . 
Similarly, we find 

\d z Sw n \ s = \A 2 5X n \ a <ec\Tr}\ s 

\ V * 5vn \* ^ c\Vs*V n \. + \A 3 5r\s + \A 3 SX n \ s 

( 3 - 49 ) , -n 

< c (\SV n z \ s + \d zz 5$ n \ s + \A 3 SX n \ s ) < ec \Tl n \ s 
=> \d z 5u n \ s <ec\K n \ s . 
Estimating the terms in 

(3.50) 

Sg n = d z {5u n ■ Vp n + u n+1 ■ V5p n ) - V 1 - ■ {5u n ■ Vv n + u n+1 ■ V6v n ) 
= 5u n z ■ Vp n + Su n ■ Vp n z + u n z +1 ■ VS P n + u n+1 ■ V5p n z 
- V x 5u n : Vv n - 5u n • VA 2 V n - V-V l+1 : V5v n - u n+1 ■ VA 2 5ip n 
using what we have computed above as 

c \Sg n \ s < e \TZ n \ s (\q\ s + \U n \ s ) + en \K n \ s {\q\ s + \U n \ s ) + e (\q\ s + \U n \ s )\K n \ s 
+ en {\q\ s + \U n \ s )\n n \ s + en \U n \ s (\q\ s + \U n \ s ) + e \U n \ s (\q\ s + \U n \ s ) 
+ en (\q\ s + \U n \ s )\K n \ s + en (\q\ s + \U n \ s )\K n \ s , 

we obtain 

(3.51) \Sg n \ s < e^c\TV\(\q\ s + \U n \ s ) < e^c \K n \ s (\q\ s + ||/||). 
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Estimating (DV n z +1 )5g n and (DTZ%)Q n in (|3~T7|) using Cauchy's formula as be- 
fore and using 13.4211 . we find 

C\K + Xn+1 < I \VXn\W n \s + | \n\s;n\G n \s 

+ \d z (6^)\ s + \d z (w"+i6v?)\ a + epK 2 \K^ n 

(3.52) < y 1^1.^1^1.(191. + 11/11) + f (\q\ s + WfWflKUn 

+ \KUn \q\ s + e\K n \(\q\ s + 11/11) + ^\n n \ s 

\K + Xn + 1 < ^S 1/4 (\q\ s + ll/H +C 6 (/i)) 2 |7^| s; „ 

The estimates for 7£™ and TV^ are identical in form, giving us 

(3.53) \n n +% ;n+1 <c 5 e l / i (\q\ s + \\f\\ + c 6 f\Tl n \ s ., n . 

So |7£ n | 5 . n is decreasing in n for s sufficiently small. Remembering that this con- 
struction is valid for n = 0, • • • , n*, |7£™| s . n is therefore smallest for n = rt» = 
rye -1 / 4 . More precisely, when 

(3.54) c 6 eV*(\q\ a + \\f\\ + c 6 ) 2 < K < 1 , 
we have 



l^ n 1.;„, <l^°l.;0 



(3-55) < |^°| s;0 exp(e- 1 / 4 7 ? logX) 

< c (\q\ s + 11/11) V 3 / 4 exp^-V^logA-) 

which is exponentially small in e (the argument of the exponential is negative since 
K < 1). In what follows we fix if, say K = 1/2, and set ?y = er/ log 2. Assuming in 
addition that e < cr/10, we can write e -3 / 4 exp(— cr/e 1 / 4 ) < c cxp(— er/e 1 / 4 ) and 

(3.56) | s .„, < c + ll/H) 2 exp(-a/ £ - 1 / 4 ). 

Remarks. A similar result for a different model was obtained in This result, 
which does not depend on Gevrey regularity (the model was inviscid), appears to 
be peculiar to the model in question and is only local in time (i.e. only up to this 
point in the present proof). 

3.3. Regularity of QGE™. We now show that the solution q(t) of 

(3.57) d t q + P < [V 4 - • (u* ■ Vv*) - d z (u* ■ Vp*)] = pA 3 q + f< 

where u* — it™* and p* = p' 1 *, is bounded independently of e for t 6 [0, To]. To 
this end we write u* = u° + u £ and p* = p° + p e where u° = (it , v°, 0), 

(3.58) v° = V ± A 3 1 q and p° = -d z A 3 1 q. 

As a preparation, we note that H3.19fa ) and (|3.23() imply, with (|3.39() . 

. ocn , |V 3 //| s <c|A 3 $*U< ce 1 / 4 (| g t + 11/11) 
(3.59) 

|V 3 ^|. < c + |A 3 X*| S + |A 3 $*| S ) < ce 1 ' 4 (\q\* + ||/||). 
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Noting that V X ■ (v° ■ Vv") - d z (v" ■ Vp°) = v° ■ Vq, we write (fTTTKTI as 

d t q + P < [v°.Vq + V 1 - ■ (u e ■ Vv°) + V 1 - ■ (v° • V« e ) + V- 1 • (u e ■ Vv e ) 



(3.60) 



d z (u" ■ Vp°) - d z (v° ■ V P £ ) - d z (u E ■ Vp e )] - fxA 3 q + /< 

2/ 



We multiply this by A^q in L (M) and estimate the resulting terms as 

Ls+i 



\(AZqy.V$\<c\q\,\v°-Vq\ a <c\q\M.-M 



(3 - 61) <^M +1+ jm-Mi-, 

and for the terms involving 

\(A s 3 q, V x ■ (u £ • V«°))| <clgl.lv- 1 • (« e • Vv°)\ s < c\q\ s k\u £ ■ Vv a 
(3.62) <c\q\ s K\u%\V 3 V \ s <c\q\ s KE 1 / i (\q\ 2 s + 

<c\€m 2 s + \\f\\), 



(3.63) \(A s 3 q,V ± -(v -Vv^))\<c\q\ 2 (\€ + 

\(A s 3 q, V 1 - ■ (rf ■ V« e ))| < c \q\ s k\ U z\ s |V 3 « e 

<ce^\q\ s (\q\ 2 s + 
and for those terms involving p, 

(3.65) \(A s 3 q,d z (u s -Vp°)\ <c\q\ 2 s (\q\ 2 s + H/II), 

(3.66) |(Agg, a z («° • Vp 6 )] < c |g| 2 + 11/11), 

(3.67) \(A s 3 q,d z (u e ■ Vf)\ < ce 1 ^ \q\ s {\q\ 2 s + ||/||) 2 . 
Putting these together, we have 

(3.68) ±\q\* + ^ | $) a +x < Cg(/x) ^(i^ + + ^ + || /|| ^ , 

where c 8 and c' are independent of e if one assumes (|3.41ll and l|3.54|l . Now let 
Q = 2 |g(0)| s + ||/||. Replacing \q\ s by Q in and l(331)l . we assume 

e < e = minjVlO, l/(16/i 2 ), 1/c 4 , ( Cl /4) 4 / 3 , l/(12 M c 3 Q) 4 , 

(3-69) l/[2c 5 (Q + ll/H +c 6 )] 4 , 1/(12 C3 ||/||) 4 / 3 , 

l/[l2c 3 c 4 (2g+||/||) 2 ] 4 ). 

It now follows from (|3.68|) that there exists a To(|q(0) | s , ||/||; /i, s, cr, M) > such 
that for t £ [0,T ], 

(3.70) |g(t)|.< 2 |g(0)| a + [[/[[. 

Remarks. It is clear from the foregoing that this result is also valid for any n < n». 
When n = (i.e. the classical quasi-geostrophic equation), one does not need the 
K cutoff to prove boundedness (cf. [3]), but this seems unavoidable for n > 0. 
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3.4. Long-time Error Bounds. Let W* = (v*,p*) = (v n *,p n *). We now seek 
to show that W — W < — W* , if initially exponentially small, remains exponentially 
small over a timescale of order one. We use the usual Gronwall-type argument. 
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Starting with the evolution equation for p < , 

d t p< - ~™< + P < (u< • Vp<) = pA 3 p< + /< 



(3.71) 

=: pA 3 p< + f< + T> , 



- P («< • Vp> + «> • Vp< + m> • Vp>) 

< _1_ f < _L T> 



we find that p = p* — p < is governed by 

(3.72) d t p ~^w + P < (u-Vp + u-Vp* +u* ■ Vp) = 5 p * + AiA 3 p + ;F>, 
where 

1 



(3.73) 5; = -ftp* + 7 «;* - P (u* ■ Vp*) + pA 3 p* + /; 



Now >S* can be expressed in terms of 7£t as follows: 

-A 3 d z S* p = d t (A 2 <$>* zz - d zz q)--A 3 d z w* + A 3 d z (u* ■ Vp* 



e 

- pA 3 d zP * - A 3 d z f< 



(3.74) 

= A 2 (D$!J - 9 22 e* + -A 3 A 2 X* + A 3 ft(u* ■ Vp*) 
- pA 3 d zP * - A 3 dJ< , 
which gives upon using the definitions (|3.14() of Q* and (|3.45() of TZ^ 

(3.75) A 3 d 2 s; = a 2 ^;. 

Similarly, we find for v, 

, d t i> + - \i) ± + Vp] + P K (u ■ Vi) + u ■ Vv* + u* ■ Vv) 

(3.76) s ' 

= MA 3 f) + S* + T> 

where 

(3.77) T> = -P < (m< • Vv> + u> ■ Vv< + u> ■ Vv>) 
and, following the computation leading to l|3.75|l . 

(3.78) S*=K, A 3 V ± ■ S* v = A 2 li; and A 3 V • S* v = A 2 TZ* X . 
Multiplying (|3~7T)| by v and $FT% by p in £ 2 (M), we find 

(3.79) ^(l*l 2 + l^)+MlV 3 «| 2 + |V 3 p| 2 ) 

= -(v,u-Vv*)-(p,u-Vp*) + (v,S* v ) + (p,S* p ) + (v^>) + (p,T> 
Noting that 

, oon , |V3«*|z-(M)<c|V30*| a <c(|?|. + |^|.)<c(|?| a + 

(3.80) . 
\V3 P *\ L ~ (M) <c\V 3 p%<c(\q\ s + \U*\ s )<c(\q\ s + 

where we have used <|3.42ll for the last inequalities, we bound the first two terms 
on the rhs of (|3.79(l by 

\{v,u- W)l < c \v\ \V 3 i,\ \V 3 v*\ Lx(M) < ^\V 3 i,\ 2 + -\vf (\q\l + ||/|| 2 ) , 

(3.81) ° £ 

\( P ,u- Vp*)\ < c \p\ \V 3 v\ \V 3 p*\ Laa(M) < ^\V 3 v\ 2 + -|p| 2 (|g| 2 + ||/|| 2 ) , 
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and use the Cauchy-Schwarz inequality for the remaining terms. This gives us 

(3.82) „1(|«| 2 + \p\ 2 ) < c(|g| s 2 + ||/|| 2 ) (|*| 2 + \p\ 2 ) + d |S*| 2 + c" \T> | 2 , 

where we have denoted |5*| 2 = |5^| 2 + |S;| 2 and \T>\ 2 = \F>\ 2 + \T>\ 2 . 

B y lTT77ra . we can replace \q(t)\ s by 2 |q(0)| s + ||/|| for t G [0,T ]. Using (1X751) 
and l|3.78[) we can then bound 

(3.83) |S*(t)| < c \TZ*\ a <c(\q(Q)\ 2 + Wff) cxp^a/e 1 / 4 ), 

also valid for t G [0, To]. 

The last term in (|3.82() can be bounded as 

l^lo < c\V 3 W<\ \V 3 W > \ < c\V 3 W\ |V 3 V^ > | 
(3-84) , ,, 

<cexp(-a/ £ 1 / 4 )|V3^|||A 3 ^|| CT . 

We now take d = 2|V 3 W | and C a = 2||A 3 W L in Lemma[T] Since |V 3 W | < 
c ||A 3 Wo|| CT by our hypothesis (i.e. a > and the rhs is finite), we have 

(3.85) \V 3 W(t)\ <2|V 3 W | and \\A 3 W(t)\\ a < 2||A 3 W L 

for t G [0,T^] where T^HAgWolU, \\f\\;a,fj,,M) = min^i,^). Thus HTm be- 
comes 

(3.86) |^>(t)| < c exp(-a/ £ 1 / 4 ) ||A 3 W || 2 CT 

fortG[0,T^]. 

Finally, writing y = \v\ 2 + \p\ 2 , we rewrite (|3.82|) as 

(3.87) < My + N 

where M is bounded independently of e, 

(3.88) M{t) = - {\q{t)\l + Wff) < ~ (l?(0)| 2 + \\ff) < -(\\A 3 q t + ||/|| 2 ) 

H ix fi 

for t G [0,T ]. As for N, $T8fy and &Ef$ give us 

N(t) = -(\S*(t)\ 2 +T>(t)\ 2 ) 

(3.89) V 

< c 9 ( M ) [(|V|Wq| + ll/ll) 2 + ||A 3 Wb||*] exp(-2 ( r/ e 1 /4) j 

valid for t G [0,T«] where T^([|ViWoL, ||/||; a, M , JVC) = min(T ,r;). 

Since W > and W are L 2 (M)-orthogonal, the hypothesis l|2.12|l implies that 

(3.90) \W Q \ 2 + \W>\ 2 < a d exp(-2a/ £ 1 / 4 ). 

Lemmanjand <|2.7[1 then imply that there is a T 3 (|| V^WoHo-, ||/||; <t, p,M) such that 
fortG[0,T 3 ], 

(3.91) \W>{t)\ 2 < 2C id exp(-2a/e 1/4 ). 

Integrating (|3.87|) and taking into account (|3.88|I - (|3.9U|) give us 

(3.92) \W{t)\ 2 < 2C id exp(-2ff/e 1 / 4 ) 

for t G [0,T 4 ] with T 4 = T 4 (|| V|Wo|| CT , \\f\\;a, ft,M). It thus follows that 

(3.93) \W>(t)\ 2 + \W(t)\ 2 < 4C id exp(-2a/ £ 1 / 4 ), 
for t £ [0, min (T 3 ,T 4 )], which is precisely l|2~T5j) . 
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